ARTICLE IN PRESS 



NP5579_proof ■ 2 September 2014 ■ 1/17 



Neuropharmacology xxx (2014) 1-17 




ELSEVIER 



Contents lists available at ScienceDirect 



Neuropharmacology 



journal homepage: www.elsevier.com/locate/neuropharm 



Neuro 

pharmacology 



Cortical EEG oscillations and network connectivity as efficacy indices 
for assessing drugs with cognition enhancing potential 

A. Ahnaou , H. Huysmans, T. Jacobs, W.H.I.M. Drinkenburg 

Dept. of Neuroscience, Janssen Research & Development, A Division ofjanssen Pharmaceutica NV, Turnhoutseweg 30, B-2340 Beerse, Belgium 




ARTICLE INFO 



ABSTRACT 



Article history: 
Received 25 April 2014 
Received in revised form 
18 August 2014 
Accepted 20 August 2014 
Available online xxx 



Keywords: 

Cognition enhancers 
EEG oscillations 
Coherent functional network 
Neurodegenerative disorders 
Animal model 
Translational biomarker 



Synchronization of electroencephalographic (EEG) oscillations represents a core mechanism for cortical 
and subcortical networks, and disturbance in neural synchrony underlies cognitive processing deficits in 
neurological and neuropsychiatric disorders. Here, we investigated the effects of cognition enhancers 
(donepezil, rivastigmine, tacrine, galantamine and memantine), which are approved for symptomatic 
treatment of dementia, on EEG oscillations and network connectivity in conscious rats chronically 
instrumented with epidural electrodes in different cortical areas. Next, EEG network indices of cognitive 
impairments with the muscarinic receptor antagonist scopolamine were modeled. Lastly, we examined 
the efficacy of cognition enhancers to normalize those aberrant oscillations. 

Cognition enhancers elicited systematic ("fingerprint") enhancement of cortical slow theta (4.5 
—6 Hz) and gamma (30.5—50 Hz) oscillations correlated with lower activity levels. Principal component 
analysis (PCA) revealed a compact cluster that corresponds to shared underlying mechanisms as 
compared to different drug classes. Functional network connectivity revealed consistent elevated 
coherent slow theta activity in parieto-occipital and between interhemispheric cortical areas. In rats 
instrumented with depth hippocampal CA1-CA3 electrodes, donepezil elicited similar oscillatory and 
coherent activities in cortico-hippocampal networks. When combined with scopolamine, the cognition 
enhancers attenuated the leftward shift in coherent slow delta activity. Such a consistent shift in EEG 
coherence into slow oscillations associated with altered slow theta and gamma oscillations may underlie 
cognitive deficits in scopolamine-treated animals, whereas enhanced coherent slow theta and gamma 
activity may be a relevant mechanism by which cognition enhancers exert their beneficial effect on 
plasticity and cognitive processes. The findings underscore that PCA and network connectivity are 
valuable tools to assess efficacy of novel therapeutic drugs with cognition enhancing potential. 

© 2014 Published by Elsevier Ltd. 



qi 1. Introduction 



A prominent property of neuronal networks is their tendency to 
engage in oscillatory rhythms. Neuronal oscillations represent a 
fundamental mechanism enabling coordinated activity during 
normal brain functioning, and are therefore an instrumental 
research target for neurological and neuropsychiatric disorders. 
Brain networks communicate with each other by means of 
neuronal oscillations at multiple temporal and spatial scales to 
integrate feature information during sensory processing and 
ensure higher order coordination of motor and cognitive functions 
(Babiloni et al., 2011; Buzsaki and Watson, 2012; Hammond et al., 
2007; Uhlhaas and Singer, 2010). Neuronal assemblies are able to 
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modulate their level of synchrony without affecting their firing rate 
(Singer, 1999), which typically occur at different frequencies both 
during sleep (Buzsaki, 2006; Steriade et al., 1993) and awake states 
(Fries et al., 2001; Pesaran et al, 2002). Oscillations in each fre- 
quency band may have a different underlying mechanism and 
subserve a different function. 

Synchronous oscillatory rhythms in the slow theta and gamma 
oscillations represent a core mechanism for sculpting temporal 
coordination of disparate brains networks during higher order 
cognitive tasks including attention and working memory (Cantero 
and Atienza, 2005; Kaiser and Lutzenberger, 2003; Lisman and 
Buzsaki, 2008; Singer, 1999). Several studies have focused on the 
origin of theta rhythm, which generators have been found in spe- 
cific hippocampal and extra-hippocampal regions (Buzsaki and 
Watson, 2012; Montgomery et al., 2008; Vertes, 2005). Oscilla- 
tions at theta frequency range are critical to establish precise 
temporal interactions to propagate and coordinate the flow of 
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Fig. 1. Time course of the effects of A/donepezil, B/rivastigmine and C/memantine on the EEG spectral patterns in fronto-parieto-occipital cortical areas (only right hemisphere is 
presented) during each 15 min block of the recording session after the administration. Mean EEG spectra were expressed as a percentage of the mean spectra in the baseline period. 
Peak changes were elicited in both slow theta (4-6.5 Hz) and gamma (32-48 Hz) frequency oscillations during 2 h after the administration. Changes from baseline are presented as 
a heat map with cold dark blue to warm red color indicates an order increase in the magnitude of oscillatory power density in particular frequency. (For interpretation of the 
references to colour in this figure legend, the reader is referred to the web version of this article.) 
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visual inspection. Only continuous waking artifact-free epochs with low-voltage fast 
EEG activity, high to moderate EMG and body activities were used in the analysis. 
Epochs with high-voltage slow cortical waves in the absence of EMG and locomotor 
activities were discarded. A notch FIR filter at 50 Hz was applied to avoid voltage 
related to power line interferences After a stable baseline recording session of 
30 min, signals from 6 brain areas were recorded for 2 h after drug administration 
(n = 8 for each condition). Continuous EEG and EMG field potentials were acquired 
at 2 kHz sample rate with an input range of ±500 mV. The signals were digitized 
with 24 bit resolution, amplified and analogue band-pass filtered between 1 and 
100 Hz. 



2.2.1. Spectral power 

Spectral density estimates were calculated using Fast Fourier Transform with 
Hanning window function; block size of 512 data points, giving 1.0 Hz resolution, 
and power was expressed as percentage of total power over 1—100 Hz. The average 
spectral density in each frequency oscillation was normalized across animals to 
obtain the full power spectrum (heat map with 1 Hz resolution) constructed from 
overlapping windows of 4 s data and plotted for baseline periods and post-drug 
periods windowing over 1—100 Hz. Frequency bands were averaged for each 
15 min session into frequency windows: Delta band (1^ Hz), Theta band (thetal: 
4-6.5 Hz, theta2: 6.5-8 Hz), Alpha band (alphal : 8-11 Hz; alpha2: 11-14 Hz), Beta 
band (betal: 14-18 Hz; beta2: 18—32 Hz), Gamma band (gammal: 32-48; gamma2: 
52-100 Hz). 

2.2.2. Principal component analysis 

The multivariate statistical principal component analysis (PCA) method was 
applied to reduce dimensionality of the datasets while retaining the variation that is 
originally present. The obtained vector component that represents the most 
important variables influencing the results diversity was used in projection tech- 
nique to look for pattern recognitions (Dimpfel, 2003). The PCA used existing 
datasets of reference drug classes obtained in the same standardized EEG experi- 
ments described in the present work. By using the fit from the PCA analysis, we can 
project results of the new experimental compound into the multi-dimensional 
component space formed by the reference compounds and we fit a cluster 
dendrogram using each compounds geometric component space mean. The 
resulting dendrogram gives an indication as to which reference compounds' class 
the newly tested compound is most similar. 

The EEG signals were measured from six brain regions in animals that received 
vehicle and drugs (low, middle and higher doses; n = 8 per condition) and spectral 



dynamics were analyzed in 9 different EEG bands during 8 time-blocks of 15 min. 
For each drug-class dataset, the doses estimated to have a therapeutic response in 
vigilance states and behavioral tests in rodent were used in the analysis. As the PCA 
requires complete observations, animals with missing values in the dataset were 
excluded from the analysis. The procedure starts with a dataset containing all 
complete vehicle observations and all complete observations at the specified dosage 
level from the reference datasets. For each animal in the dataset, there were 432 
variables, representing every possible combination between time-point, brain re- 
gion and EEG bands. The basic effect of the PCA projection is to project all the var- 
iables in the dataset to a small set of components that account for the greatest 
proportion of the variance in the data. Using these 432 variables, a principal 
component analysis was fitted to all observations to the specified number of com- 
ponents. Neither the drug name nor the drug classes were passed to the PCA 
function. The component space is made purely as a function of the EEG information 
contained in each animal's 432 variables. The drug class is strictly used for coloring 
the PCA projections and dendrogram. Also note that the test compound dataset was 
not used to form the PCA projections and projections for the reference compounds 
were identical across runs of the PCA application. 

The variable loadings from the PCA projections were saved and applied to the 
newly test compound dataset. Then, the drug's geometric center was computed by 
averaging the locations of all observations (animals) for that drug. This averaging 
will have dimensionality equal to the number of specified principal components. 
Thus, for a PCA projection fit that is using 4 principal components, the geometric 
center will be a set of coordinates in 4-dimensional space. These components were 
fit into hierarchical clustering dendrogram using the Ward's linkage, which helps to 
see and understand the group structures in the data. Here, we have used as our 
unknown compound "Rivastigmine" dataset obtained in a standard EEG experiment, 
which serves as a good example of the application's abilities. 

2.2.3. Functional network activity 

Functional network activity between cortical regions was computed by coher- 
ence, which is commonly used to probe the integrity of cortical neural pathways and 
index the functional coupling between different cortical structures at various fre- 
quencies (Thatcher, 2012). The amplitude of EEG coherence ranges between 0 and 1: 
a low value indicates no similarity between the two signals, whereas values close to 
1 indicate a high similarity between two time series signals. Coherence between the 
time series in different electrodes was quantified at different time points before and 
after administration of drugs using a generalized additive mixed model where the 
coherence over the frequencies was fitted using a cubic spline for the different 
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treatment groups and the heterogeneity between the animals (Priestley, 1981; 
Q3 Chatfield, 1975, 2000). Subsequently, coherence between each pair of electrodes 
were concatenated every 4 s and aggregated over 15 min in each frequency bands for 
each animal, time point, and treatment group. Particular frequency bands were 
selected for comparison between groups based on the most prominent oscillatory 
and coherent activities across different brain regions. To reduce the impact of vol- 
ume conduction of signals propagated from common generators, computation of all 
edges in which the maximum absolute value of coherence occurred at zero time lag 
between time series were removed. 

2.2.4. Epidural versus depth EEC recordings 

A separate group of 20 rats was stereotaxically instrumented with intracortical 
(Ml : AP +2 mm, L ± 2 mm, V +2 mm from the Bregma) as well as intra-hippocampal 
electrodes (CA1: AP -3.14 mm, L ± 1.8 mm, V +2.8 mm; CA3: AP -4.5 mm, 
L ± 3.8 mm, V +4 mm from Bregma, respectively) in order to compare functional 
coherent activity between cortical (epidural) and hippocampal (depth) networks. 
Donepezil at the dose of 3 mg/kg shown in earlier studies to elicit peak coherent 
activity in epidural recordings was used in this study. 

2.3. Drugs 

All drugs were synthesized in Janssen Research & Development laboratories. 
Donepezil (0.3, 1.0 and 3.0 mg/kg), rivastigmine (0.3, 1.0 and 3.0 mg/kg), tacrine 
(0.16, 2.5, 10 mg/kg) galantamine (0.16, 2.5, 10 mg/kg), memantine (1, 3 and 10 mg/ 
kg), scopolamine (0.16, 0.64 and 2.5 mg/kg) and amphetamine (0.16, 0.64, 2.5 mg/kg) 
were formulated in sterile saline solution, and were administered subcutaneously 
(s.c.) in a volume of 10 ml/kg body weight. Saline was administered in control 
groups. All drug classes used in the database were used at the doses that showed 
therapeutic efficacy in neurophysiological and behavioral studies in rodents. The 
database includes antidepressants (desipramine, nortryptiline, paroxetine, fluoxe- 
tine, trazodone buproprion); Anxiolytics (ipsapirone, oxazepam, buspirone, Zolpi- 
dem); psychostimulants (amphetamine, PCP, ketamine); antipsychotics (olanzapine, 
risperidone, clozapine, paliperidone, haloperidol); mood stabilizers (lamotrigine, 
valproate). 

2.4. Statistics 

In the dose response experiments, values of consecutive 4-s epochs were 
averaged for a 30 min period (baseline) followed by 15 min periods (treatment) 
averages for the remaining time after the administration of vehicle or test com- 
pound. In the pharmacological reversal challenge experiments, EEG oscillations 
were averaged for 30 min before the muscarinic cholinergic antagonist (baseline), 
and then for every 15 min after combined treatment with cognition enhancers. 
Drug-induced changes in EEG spectral power were calculated in blocks of 15 min for 
a total of 2 h as the ratio of mean spectral power obtained following the adminis- 
tration of test drug versus the mean spectral power obtained during 30 min baseline 
period. This procedure allows for assessment of drug-induced changes in different 
cortical networks and EEG power expressed in each frequency band as a percentage 
of the original power, compared to the vehicle condition. Time course EEG spectral 
changes were submitted to a two-way multivariate ANOVA for repeated measures 
with two main factors (treatment and period of recordings), followed by a one-way 
ANOVA with main factor treatment for each of the 15 min period. Data were pre- 
sented as the means ± S.E.M. The strength of association between locomotor activity 
and EEG slow theta oscillation (4—6.5 Hz) was assessed with linear Spearman cor- 
relation analysis. The aggregated coherences as well as the changes from baseline 
were analyzed per frequency band using an ANOVA with time, group and its 
interaction as covariates while taking into account the heterogeneity between the 
different animals (Verbeke and Molenberghs, 2000). P-values from the associated F- 
tests were reported as well as Least squares means from the corresponding ANOVA 
models. 



3. Results 

3.1. Cognition enhancers enhanced EEG slow theta (4—6.5 Hz) and 
gamma (32—48 Hz) oscillations 

We present the effects of different cognition enhancers on EEG 
oscillations. AChEIs donepezil, rivastigmine, tacrine, galantamine 
and the NMDA antagonist memantine elicited common changes in 
EEG network oscillations across different brain areas. A concomi- 
tant dose-dependent increase in the slow theta (4—6.5 Hz) and 
gamma (32—48 Hz) oscillatory patterns were found in different 
brain regions (Fig. 1 A— C heat maps of the right side hemisphere are 
displayed for donepezil and rivastigmine and memantine). 

Quantitative assessment of slow theta (4—6.5 Hz) and gamma 
(32—48 Hz) frequency power for all cognition enhancers is 



indicated in Fig. 2A and B. A much pronounced effect on slow theta 
activity (4—6.5 Hz) was observed particularly in parietal and oc- 
cipital cortical areas, whereas an enhanced effect on slow gamma 
(32—48 Hz) occurred in all cortical areas. However, no major effect 
was found in the higher gamma frequency range (52—100 Hz). 
Thus, oscillatory activities in the slow theta and gamma frequency 
range appear to be particularly responsive in the mechanism of 
action of cognition enhancers. 

3.2. Correlation between slow theta oscillations (4—6.5 Hz) and 
locomotor activity 

It was assessed whether the time course of the increase in EEG 
slow theta activity (4—6.5 Hz) oscillations were associated with 
changes in motor behavior following the administration of cogni- 
tion enhancers. Analysis of motion activity at different time points 
showed that administration of all cognition enhancers' was asso- 
ciated with a reduction in activity levels (Fig. 2C) as opposed to the 
effects of the classical psychostimulant amphetamine. 

The results of a Spearman rank correlations analysis at 30 min 
after the administration of cognition enhancers between the slow 
theta (4—6.5 Hz) and activity levels are shown in Fig. 2D, upper 
panels. The correlation analysis showed that slow theta (4—6.5 Hz) 
was negatively correlated with activity levels (donepezil: r = -0.85, 
p < 0.001; rivastigmine: r = -0.57, p < 0.001; tacrine: r = -0.63, 
p < 0.001; galantamine: r = -0.66; memantine: r = -0.43, 
p < 0.05). The scatter plot in Fig. 2D, bottom panels shows a highly 
significant negative association between EEG slow theta activity 
(4—6.5 Hz) and activity for each time intervals across all treatments 
(Spearman rank correlation analysis for all time points: rS ranges 
between -0.35 and -0.97, all P-values < 0.05). 

However, the Spearman rank correlations analysis at 30 min 
after amphetamine administration showed a positive correlation 
between slow alpha activity (8-11 Hz) (r = 0.83, p < 0.001) and 
motor activity, which was maintained across all time intervals of 
the recording session (rS ranges between 0.7 and 0.84, all P- 
values < 0.001). 

3.3. Principal component analysis 

Principal component analysis was used to identify common 
compact clusters that correspond to shared underlying mecha- 
nisms ("fingerprinting'). An unknown compound "Rivastigmine" 
dataset (black filled circles; Fig. 3) was used, which was not 
considered in the original learning dataset, and was projected 
into the space defined by our historical drug-classes datasets. 
Here, the first two components account for 20.8% and 17.1% of the 
variance in the dataset, respectively (Fig. 3, left panel). Using the 
four components, we fit a hierarchical clustering and visualize 
the clustering with a dendrogram. The results show that riva- 
stigmine was grouped in the common spatial pattern of cognitive 
enhancing drugs, indicating that its effect on EEG oscillations 
over time was most closely matched with the cognition en- 
hancers donepezil, tacrine, galantamine, memantine (Fig. 3, right 
panel). 

3.4. Parallel coherent patterns in epidural and depth cortical 
structures 

To examine whether frontal, parietal and occipital networks 
were functionally related, we computed coherence indices in 
fronto-parieto-occipital cortical areas following the administration 
of donepezil at different doses. This revealed consistent increases in 
slow theta (4—6.5 Hz) and gamma (32—48 Hz) activities accom- 
panied by coherence increase in those frequency bands in 
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Fig. 2. Changes in the power of A/slow theta (4-6.5 Hz), B/slow gamma (32-48 Hz) frequencies during the first 2 h after the administration of donepezil, rivastigmine, tacrine, 
galantamine and meinantine (only right hemisphere is displayed). The red dots on the brain diagram mark the brain region in which the drug effects are displayed. Color coded bars 
underneath the curves indicate intervals in which oscillatory activity difference differed from vehicle Mixed model ANOVA. C/Activity levels during the first 2 h after the 
administration of donepezil, rivastigmine, tacrine and galantamine. Note that the classical psychostimulant amphetamine was included to display differences in motion levels. D/ 
Spearman correlations (upper panels) at 30 min after the administration of different cognition enhancers between the slow theta oscillations (4-6.5 Hz) and locomotor activity, and 
(lower panels) R-Values plotted for each time block of 15 min. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this 
article.) 
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Fig. 2. (continued). 
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Fig. 3. A/A scatter plot of the principal component analysis for all drug classes tested. The PCA analysis was based on 432 variables (six brain regions times 8 for 9 frequency bands). 
The space projection distance biplot displayed the mean adjusted observations of the main two obtained principal components. The colored symbols in the figure legend 
correspond to the six drug classes and vehicles. Different drugs classes are coded by different colors. B/Dendrogram drug clusters. The level at which branches merge (relative to the 
"root" of the tree) is related to their highest similarity: the height of the branch points indicates how similar or different are from each other: the greater the height, the greater the 
difference between clusters. Note that the drug Rivastigmine, which was not part of the initial learning datasets showed high similarity with other cognition enhancers in the left 
part of the plot. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) 



distributed fronto-parieto-occipital networks (Fig. 4A, coherent 
pairs left panels). 

The presence of slow theta (4—6.5 Hz) and gamma (32—48 Hz) 
frequency oscillations suggests that cortical networks might 
spontaneously generate those oscillations or be derived from vol- 
ume conduction of the signals propagated from common source. 
Therefore depth EEG recordings within the motor cortex, CA1 and 
CA3 after the administration of vehicle (n = 10) and of donepezil at 
the dose of 3 mg/kg (n = 10) was performed. In accordance with 
epidural recordings, the presence of prominent slow theta 
(4—6.5 Hz) and gamma (32—48 Hz) power associated with peak 
phase coherent activity was found. Power spectra derived from 
intracortical and intrahippocampal sites in rats implanted with 
depth electrodes displayed peak coherent slow theta (4—6.5 Hz) 
and gamma (32—48 Hz) oscillatory activities. Cortical motor cortex 
M1/IV12 areas displayed marked gamma oscillations in all leads 
whereas CA1 and CA3 exhibited slower theta (4—6.5 Hz) rhythm. 
Donepezil (3 mg/kg) elicited strong cortico-hippocampal slow 
theta (4—6.5 Hz) and gamma (32—48 Hz) coherent activities (Fig. 4, 
coherent pairs right panels). This pattern of synchronized activity of 
the hippocampus and cortical structures parallels the pattern seen 
in rats instrumented with epidural electrodes. These data suggest 
those epidural network oscillatory activities in the slow theta 
(4—6.5 Hz) and gamma (32—48 Hz) range are unlikely derived from 
the same source as removal of the edges with maximal absolute 
coherence values occurring at zero time lag excluded spurious 
edges attributable to volume conduction. 

Analysis of functional connectivity changes as evolving over 
time, showed maintenance of stable coherent activity in slow theta 
(4—6.5 Hz) frequency oscillations with consistent contralateral 



coherent activities (Fig. 4B). In addition, a strong dose-dependent 
decrease in coherent slow alpha (8—11 Hz) was observed in 
fronto-parieto-occipital pairs. 

3.5. Effects of scopolamine on EEG oscillations and network 
connectivity 

To model EEG correlates of cognition deficits observed in brain 
disorders involving cognitive functions, we dose-dependently 
influenced the cholinergic transmission by the muscarinic cholin- 
ergic receptor antagonist scopolamine (0.16, 0.64 and 2.5 mg/kg). 

Scopolamine evoked a dose-dependent leftward shift towards 
slow EEG network oscillations. A peak of functional effect 
expressed most clearly in the increased slow delta activity in 
different brain areas was found at about 15 min after administra- 
tion of the middle and highest doses and lasted till at least 2 h post- 
treatment (Fig. 5). A striking and robust network coherent delta 
activity was found in fronto-parieto-occipital sites, whereas 
coherence in fast oscillatory rhythms was strikingly reduced (data 
not shown). 

3.6. Efficacy of the cognition enhancers to attenuate scopolamine- 
induced aberrant EEC oscillations and network connectivity 

Quantitative evaluation of the efficacy of cognition enhancers to 
reverse abnormalities in EEG network connectivity associated with 
disorganized behavior was performed in scopolamine-treated rats. 
The cognition enhancers tested in this study were potent to 
significantly reverse the leftward shift and increased coherent slow 
EEG delta oscillation synchrony in scopolamine-treated rats, in 
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Fig. 4. A/Network coherence changes after subcutaneous administration of vehicle and donepezil in epidural recording (left panels; n = 8 for each condition) and intracortical 
recording (right panels; n = 10 for each condition). Donepezil at 3 mg/kg elicited similar synchronized slow theta (4-6.5 Hz) and gamma (32-48 Hz) pattern in epidural and 
intracortical recording in comparison with that of the control group. B/Functional network coherence changes from baseline in each frequency band of interest derived from 
consecutive 4s epochs concatenated and averaged in 15 min blocks. The width of the edges between pair electrodes is drawn proportional to the weight of changes in coherent 
activity from baseline (0-50%, 50-100% and 100-150%). Red color represents increases and blue color indicates decreases. Average networks derived from the first session 
recording session in each pharmacological condition and for each frequency band of interest 6 (0.5-4 Hz); 91 (4.5-6 Hz); 62 (6.5-8 Hz); a\ (8.5-11 ); <x2 (11.5-13); pi (13.5-18 Hz); 
P2 (18.5-30 Hz); yl (30.5-50 Hz); y2 (50.1-70 Hz), y3 (70.1-100 Hz). Time-evolving dominant coherent activity in the frontal-parietal and occipital structures indicated most 
persistent edges scaled to 01 frequency during waking. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) 
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Fig. 5. A/Time course of the effects of subcutaneous administration of scopolamine (0.16, 0.64, and 2.5 mg/kg) on the EEG spectral patterns in fronto-parieto-occipital cortical areas 
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location dependent manner (Fig. 6A, only Donepezil displayed). The 
findings suggest restoration of abnormalities in the network coor- 
dination to its large-scale synchronization form. 

Although no behavioral cognitive task was considered in these 
studies, the leftward shift in EEG neuronal oscillation likely indi- 
cated a depression in brain activity, which could correlate with 
impairment in cognitive processing. 

4. Discussion 

The present studies investigated the effects of cognition en- 
hancers, which are approved for symptomatic treatment of patients 
with mild-moderate AD, on EEG oscillations in freely behaving rats 
and their potency to attenuate aberrant oscillations and network 
connectivity in the scopolamine-induced cognitive deficits model. 
Three major results are reported: i) EEG field potentials following 
cognition enhancers manifested common changes, notably a 
concomitant enhancement of oscillations and connectivity in slow 
theta (4—6.5 Hz) and gamma (32—48 Hz) in different brain net- 
works associated with reduced motion activity, ii) Scopolamine 
promoted a consistent leftward shift coherent activity confined to 
the delta-range, iii) pre-treatment with cognition enhancers as 
shown with donepezil attenuated the leftward shift in network 



connectivity dominated by enhanced coherence in slow EEG 
oscillations. 

4.1. Translational value of the EEC recording sites 

The rat's frontal cortex and specifically the medial frontal cortex, 
is usually considered to be homologous to the dorsolateral pre- 
frontal cortex of primates, which consists of cortex homologous to 
human anterior cingulate cortex (ACC) and dorsolateral prefrontal 
cortex (Seamans et al., 2008; Wallis, 2012). The cortex is hierar- 
chically organized and some cortical areas serve higher-order 
integrative functions that are neither purely sensory nor purely 
moto. These higher-order areas of the cortex, known as association 
areas, serve to associate sensory inputs to motor response and 
perform those mental processes that intervene between sensory 
inputs and motor outputs. The cortex of each hemisphere contains 
one primary motor area (Ml) and primary somatosensory areas. 
The Ml is not a simple transmission relay between the non- 
primary motor areas (that anticipate, prepare) and the spinal 
cord, but rather a real crossroad playing an important role in 
cognitive motor integration as it could be active without any motor 
tasks (Carlsson et al., 2000; Saper et al., 2000). The primary so- 
matosensory cortex and parietal association cortex provide infor- 
mation about the ongoing motor response. Because the multimodal 
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Fig. 6. A/Time course of the effects of combined subcutaneous administration of vehicle/vehicle, vehicle/scopolamine (0.16 mg/kg), donepezil (3 mg/kg)/vehicle and of donepezil 
(3 mg/kg)/scopolamine (0.16 mg/kg) on the EEG spectral patterns in fronto-parieto-occipital cortical areas (only right hemisphere is presented) during each 15 min block of the 
recording session. Mean EEG spectra were expressed as a percentage of the mean spectra in the baseline period. Peak changes were elicited in both slow delta (0.5-4 Hz) frequency 
oscillations during 2 h after the administration. Changes from baseline are presented as a heat map with cold dark blue to warm red color indicates an order increase in the 
magnitude of oscillatory power density in particular frequency. B/Changes in peak delta frequency during 2 h after the combined treatments (only right hemisphere is displayed, 
n = 8 for each condition). Color coded bars underneath the curves indicate intervals in which oscillatory activity difference differed from vehicle. Mixed model ANOVA. C/Time- 
evolving dominant coherent activity in the frontal-parietal and occipital structures indicated the potential of donepezil to attenuate scopolamine-induced most persistent 6 fre- 
quency edges. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.) 
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1 association areas integrate sensory modalities and link sensory cognition enhancers appear to influence cortical activation in these 66 

2 information to the planning of movement, they are thought to be regions at lower motor activity. 67 

3 the anatomical substrates of the highest brain function in human's Although the theta oscillatory rhythm was less pronounced in 68 

4 conscious thought, perception, and goal-directed action. Hence, the frontal sites as compared to central and posterior recording sites, 69 

5 Ml and primary somatosensory cortical recording sites were there were increased stable functional slow theta (4—6.5 Hz) 70 

6 explicitly chosen in the present work. coherent edges in bilateral frontal, parietal and occipital regions. 71 

7 Spectra based coherence derived from depth recording signals in 72 

8 4.2. Relevance of theta and gamma oscillations and connectivity in donepezil treated animals exhibited a coherent pattern, which 73 

9 the mechanism of action of cognition enhancers closely resembles the patterns of synchronized activity prevalent 74 

10 at epidural recording sites. Thus, it is presumable that the 75 

11 All cognition enhancers consistently enhanced cortical oscilla- increased theta coherence in fronto-parieto-occipital cortical re- 76 

12 tions and coherent networks activity in the slow theta (4—6.5 Hz) gions was due to temporal coupling of local networks activities 77 

13 associated with lower motor activity. generated in the underlying cortical structures i.e. the hippo- 78 

14 Dynamic neural oscillations are fundamental mechanisms for campus and cortex to facilitate attention and encoding process. On 79 

15 flexible communication between brain circuits to integrate sen- a related note, it is worth mentioning that the hippocampus in 80 

16 sory information and to coordinate higher-order cognitive and rodents is located under the cerebral cortex, whereas in primates 81 

17 motor functions, suggesting that these coordinated patterns of it is located in the medial temporal lobe, underneath the cortical 82 

18 activity could be a crucial target in cognition research (Varela surface containing two main interlocking parts such as the 83 

19 et al., 2001; Womelsdorf et al., 2007). Neuronal oscillations in Amnion's horn and the dentate gyrus. Notwithstanding this 84 

20 the theta and gamma range can be measured at different anatomical disparity, functional neuroimaging and large scale 85 

21 anatomical scales, ranging from invasive recordings of small microelectrode recording cognitive or sleep studies support the 86 

22 neuronal populations till epidural regional networks. Synchroni- organization of networks, which provides an optimal functional 87 

23 zation of networks at theta rhythm establishes precise temporal brain state for effective hippocampus-to-neocortex processing and 88 

24 interactions between disparate brain areas, whereas the gamma transfer of information (Burianova et al., 2012; Cohen, 2011; 89 

25 rhythm establishes synchronization with great precision in local Ferrara et al., 2012; Le Van Quyen et al., 2010; Staresina et al., 90 

26 cortical networks. 2013; Wagner et al., 2010). 91 

27 Synchrony in the theta frequency oscillations represents one of All cognition enhancers elicited cortical oscillations and 92 

28 the best-studied neuronal activities in mammalian's hippocampus, coherent activity in the slow gamma (32—48 Hz) frequency across 93 

29 which is implicated in the top-down control of cognitive functions different cortical areas. Cortical gamma oscillations have received a 94 

30 (Vertes, 2005). Theta oscillations, which dominates the EEC great interest because of their role in sensory and cognitive func- 95 

31 recorded from the hippocampal formation and the medial septum, tions such as selective attention (Doesburg et al., 2008; Fries et al., 96 

32 acts as a pacemaker for the generation of this neuronal rhythm, 2001; Lakatos et al., 2008; Womelsdorf and Fries, 2006; 97 

33 which drives prefrontal cortical networks into theta oscillatory Womelsdorf et al., 2007), short-and long-term memory 98 

34 rhythm (Thierry et al., 2000). The strong theta rhythm synchronizes (Herrmann et al., 2004; Pesaran et al., 2002; Wu et al., 2008) and 99 

35 medial prefrontal cortex neurons and strengthens synaptic links multisensory integration (Lakatos et al., 2007). The primary gen- 100 

36 between spatially distributed cortical areas, and therefore facili- erators of gamma oscillations include the cortex, where subsets of 101 

37 tates the transfer of hippocampal information to neocortex during inhibitory GABAergic interneuron circuits that modulate gluta- 102 

38 learning and memory (Buszaki, 2002; Paz et al., 2008). It has been matergic pyramidal cell activity, are thought to play a primary role 103 

39 repeatedly demonstrated that perceptive, cognitive, and motor (Gray and McCormick, 1996; Traub et al., 1999). However, gamma 104 

40 processes are associated with enhanced EEC coherence in the theta activity is also prevalent in the hippocampus formation, where it is 105 

41 frequency oscillations between prefrontal and posterior association postulated to assist cognitive functions such as encoding and 106 

42 cortex and hippocampus (Molle et al., 2004; Seemuller et al., 2012; retrieval of memory traces (Bragin et al., 1995; Chrobak and 107 

43 Schmidt et al., 2013; Womelsdorf and Fries, 2006). Coherence is Buzsaki, 1998). The hippocampal gamma activity is segregated 108 

44 commonly used to probe the integrity of cortical neural pathways into two distinct frequency rhythms, which are believed to rise 109 

45 and index the functional coupling between different cortical from 2 distinct sources: in the CA1 nuclei, the slow gamma ranging 110 

46 structures at various frequencies (Thatcher, 2012). In the present from 25 to 50 Hz and the fast one from 65 to 140 Hz. The fast 111 

47 study, all cognition enhancers modulated the slow theta (4—6.5 Hz) gamma oscillations are entrained by the entorhinal cortex, which 112 

48 oscillations and its coherent activity, which may be important for couple CA1 to inputs from the medial entorhinal cortex, an area 113 

49 memory and cognitive process. The strong and stable coherent that provides information about the animal's position. In contrast, 114 

50 slow theta (4—6.5 Hz) frequency between cortical areas may be the slow gamma oscillations are generated within the hippocampal 115 

51 associated with enhanced functional connectivity between frontal network and couple CA1 to CA3, a hippocampal subfield essential 116 

52 cortex and hippocampus for gating information flow and cortical for storage of information input from entorhinal cortex (Bragin 117 

53 processing. et al., 1995; Colgin et al., 2009; Csicsvari et al., 2003). 118 

54 Changes in theta activity are strongly associated with changes in Modulation of the slow and fast gamma frequency oscillations 119 

55 motor behavior in rats (Oddie and Bland, 1998; Vanderwolf, 1969). plays differential roles during specific cognitive processes. For 120 

56 To rule out the possibility that EEC changes reflect the central ef- instance, during the encoding phase of the delayed-match-to- 121 

57 fects of cognition enhancers rather than changes in motor behavior, sample cognitive task, Vidal et al. (2006) examined the influence 122 

58 we investigated the changes in motion activity and whether the of attention-related information processes (subset of item differed 123 

59 time course of these alterations were correlated. In contrast to in color: attention) and visual grouping of the display (homoge- 124 

60 classical psychostimulants' effect on locomotor activity levels, the nous: grouping condition) on gamma oscillations. The visuo- 125 

61 cognition enhancers induced slow theta (4—6.5 Hz) when motor perceptual processes during attention focusing properties of the 126 

62 activity is decreased revealing dissociation between theta power stimulus required slow gamma frequency oscillations, whereas 127 

63 and motor activity. In addition, multiple spearman correlations perceptual grouping of the stimulus modulated the fast gamma 128 

64 analysis confirmed a strong negative association of slow theta oscillations, which argues for functional specialization of both 129 

65 (4—6.5 Hz) with low activity levels. These findings suggest that oscillatory rhythms in cognitive processing. 130 
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1 Our in vivo findings indicate that cognition enhancers elicited placed adjacent to the reference cognition enhancers. It is worth 66 

2 slow gamma (32—48 Hz) oscillations, which may facilitate the noting that in each experiment which contains a vehicle and three 67 

3 interaction between short distance brain networks for attentional doses (low, middle and high), there was no standardized dosing 68 

4 and information storage process. schemes between experiments and these dosages were not recor- 69 

5 It is important to note that all cognition enhancers elicited large ded in the data (e.g. 15 mg/kg), nor were standardized to an EC50 or 70 

6 regional differences in low theta oscillations being more pro- IC50 values. Thus, similar drugs might not be classified together 71 

7 nounced in the central and posterior recording sites as compared to because of a difference in the dosing scheme used. Secondly, the 72 

8 the frontal cortex (Fig. 2A). Simultaneous multisite recordings from data do not provide indication of the Tmax for each drug. Since each 73 

9 the hippocampus and the frontal cortex in behaving rats and mice of the variables in the PCA projection are simply independent 74 

10 demonstrated that hippocampal theta oscillations control the firing combinations of time-point, brain region and EEC frequency band, 75 

11 of prefrontal neurons, thereby delivering the temporal coordina- it is possible that two similar compounds would not be clustered 76 

12 tion of oscillating networks accounting for information transfer and together due to a difference in their difference in Tmax. Never- 77 

13 storage (Sirota et al., 2008; Wierzynski et al., 2009). The activating theless, the PCA identified compact clusters on the two main var- 78 

14 cholinergic and glutamatergic drive directly interact with this iables with common spatial patterns for cognition enhancing drugs 79 

15 oscillatory activity by modulating the hippocampal theta rhythm displayed. The results showed the utility of PCA to keep funda- 80 

16 and shifting the frontal cortical activity from slow to fast oscilla- mental features of variables to summarize data into compact 81 

17 tions (Metherate et al., 1992; Manns et al., 2000, 2003). Here, all clusters that correspond to common underlying mechanisms, and 82 

18 cognition enhancers increased low gamma oscillations in frontal suggest that PCA technique holds a great capacity to successfully 83 

19 cortex areas (Fig. 3A). This observation suggests that hippocampal reveal oscillatory pattern classification. 84 

20 theta entrainment of distant frontal cortex neurons into low 85 

21 gamma rhythm represent dynamic functional binding mechanisms 4.4. Scopolamine-induced slow EEG network oscillations 86 

22 of widely distributed cortical assemblies, essential for cognitive 87 

23 processing. Scopolamine elicited a consistent leftward shift in the network 88 

24 In addition, donepezil decreased coherent slow alpha (8—11 oscillatory pattern, which was often accompanied by persistent 89 

25 activity) in fronto-parieto-occipital recording pairs. increased coherent activity in slow frequency oscillations. 90 

26 The thalamo-cortical circuits generate a large network oscilla- Brain cholinergic neurons are believed to be involved in mem- 91 

27 tory activity in the alpha frequency range during cortical opera- ory and cognitive functions, and extensive research in both animal 92 

28 tions. However, this rhythm may also emerge during loss of and human demonstrated a key causal role of disruption of the 93 

29 attention in immobile state or during failure to process sensory cholinergic system in memory and cognitive decline (Klinkenberg 94 

30 input, and may serve to functionally disengage and reduce the and Blokland, 2012). Memory and cognitive functions are known to 95 

31 processing capabilities of a given brain region. The alpha activity is decline with advancing age, in mild cognitive impairment and AD. 96 

32 decreased in engaged brain region whereas it increased in disen- A major brain area of neuronal cell loss is the cholinergic nucleus 97 

33 gaged regions (Klimesch et al., 2007). In general, theta and gamma basalis of Meynert (Whitehouse et al., 1982), which sends major 98 

34 activities are accompanied by reduced synchrony in network alpha cholinergic input to the hippocampus and cortex, both areas of 99 

35 rhythm. Although the controversial views on translatability of which are highly associated with attention and memory formation 100 

36 alpha oscillatory rhythm between human and rodent, the present and storage. Imaging studies revealed a consistent decreased ac- 101 

37 results suggest that the temporal enhancement of oscillations and tivity in medial temporal lobe structures such as parahippocampal 102 

38 coherent activity in theta and gamma frequencies induced by and hippocampal regions (Celone et al., 2006; Machulda et al., 103 

39 donepezil may damp cortical alpha activity and engage those 2003; Bradley et al., 2002; Matsuda, 2001). Diffused and slowed 104 

40 cortical regions. EEG with peak in slow delta/theta network activity, higher coher- 105 

41 However, a common limitation to epidural EEG based connec- ence in delta and reduction in fast frequency oscillations reflect the 106 

42 tivity analysis may be biased by volume conduction of signals major changes in cortical activity seen in patients with senile de- 107 

43 propagated from a common generators. The model of computation mentia or AD (Babiloni et al., 2009, 2011; Gianotti et al., 2008; 108 

44 and removal of edges with maximal absolute coherence values Huang et al., 2000; Jeong, 2004; Lehmann et al., 2007; 109 

45 occurring at zero time lag described in human recordings (Chu Ponomareva et al., 2003; Rossini et al., 2006; Yener and Basar, 110 

46 et al., 2012) and used here, likely excluded spurious edges attrib- 2013). This loss of cholinergic drive has formed the basis of the 111 

47 utable to volume conduction suggesting that epidural network cholinergic hypothesis of AD, which provided the use of scopol- 112 

48 connectivity in the slow theta (4—6.5 Hz) and gamma (32—48 Hz) amine as a model of cognitive deficit. Scopolamine is a non- 113 

49 range are unlikely derived from the same source. These results are selective muscarinic receptor antagonist, which blocks the stimu- 114 

50 particularly interesting and could prove that functional network lation of post-synaptic receptors that may consequently impair the 115 

51 connectivity measures in epidural recordings may provide sensitive activity of hippocampal and cortical fast network oscillations 116 

52 indices of functional interregional network dynamics. leading to memory deficits associated with enhanced slow EEG 117 

53 oscillations in delta frequency bands (Ebert and Kirch, 1998; 118 

54 43. PCA revealed compact cognition enhancers cluster Chauviere et al., 2009; Dimpfel, 2005; Kikuchi et al., 1999). How- 119 

55 ever, the scopolamine model presents some limitations as it fails to 120 

56 The principal component analysis captured a common clus- recreate the cardinal pathological aspects and the progressive 121 

57 tering for cognitive enhancing as compared to different drug degenerative nature of AD (Van Dam and De Deyn, 2006). Despite 122 

58 classes. this caveat, the scopolamine model remains a valuable pharmaco- 123 

59 In the present work, it was noted that the test dataset (riva- logical screening tool, which offers opportunities in clinical and 124 

60 stigmine) served as a good example to assess the application's translational research for discovery and characterization of poten- 125 

61 abilities and to project exactly to where the reference drug lies in tial cognition-enhancing drugs. 126 

62 the principal component space. When looking at the distances The present results in rats demonstrate that acute disruption of 127 

63 between our test dataset and the reference drugs, rivastigmine was cholinergic transmission transiently produced abnormalities in 128 

64 the first closest drug to acetylcholinesterase inhibitors and the cortical network activities similar to those described in neurode- 129 

65 cluster dendrogram indicated that the test drug was correctly generative disorders i.e. decreased cortical oscillations and 130 
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1 networks connectivity in high frequency oscillations, and shifted to modulation of regional and local circuit's muscarinic and gluta- 66 

2 connectivity towards slow frequency rhythm. The findings that matergic receptors. 67 

3 scopolamine-induced disruption of cortical and subcortical oscil- In general, theta rhythmicity organizes gamma frequency dy- 68 

4 lations and networks connectivity further support the pivotal role namics of neuronal populations (Bragin et al., 1995; Colgin et al., 69 

5 of acetylcholine in the coordination of distributed network activ- 2009; Mizuseki et al., 2009; Sirota et al., 2008). The increased 70 

6 ities. The reduced coherence in fast oscillations associated with the slow gamma oscillations during the slow theta-entrained activity 71 

7 shift into peak coherent activity in delta frequency range in may constitute an important aspect of the communication among 72 

8 response to scopolamine reflects a cortical component of cognitive brain networks to facilitate the coupling of cortical-hippocampal 73 

9 impairment. CA3-CA1 networks in order to convey information in the memory 74 

10 stores. These experimental results and theoretical considerations 75 

11 4.5. Cognition enhancers attenuated scopolamine-induced leftward suggest the relevance of slow theta and gamma network oscilla- 76 

12 shift of EEC oscillations and peak slow coherent activity tions and modulation of these two oscillatory rhythms may, in part, 77 

13 underlie the pro-cognitive actions of the cognition enhancers in the 78 

14 The cognition enhancers effectively attenuated the early stages on AD where sufficient numbers of cholinergic and 79 

15 scopolamine-induced aberrant leftward oscillatory shift accompa- glutamatergic projection neurons are still available. 80 

16 nied with strong coherence in slow delta activity in different brain Overall, the present findings are encouraging as they suggest 81 

17 sites. that changes in cortical oscillations and networks connectivity can 82 

18 Acetylcholine system is known to be crucial for triggering theta be relevant in the mechanism of action of cognition enhancers and 83 

19 and gamma oscillation in hippocampal slices (Fellous and the preservation of cognitive processes. Functional network con- 84 

20 Sejnowski, 2000; Traub et al., 2000) and interaction of both nectivity measures in slow EEG delta as well as lower theta and 85 

21 rhythms is weakened following muscarinic receptor blockade in gamma oscillations in conscious animals may prove to be useful 86 

22 mice (Hentschke et al., 2007). Preclinical and clinical experiments efficacy indices for translational research, and provide a major 87 

23 were designed to examine the utility of the scopolamine-induced opportunity for drug discovery of drug candidates with potential 88 

24 cognitive impairment model in predicting pharmacodynamics pro-cognitive efficacy. 89 

25 signals of putatively pro-cognitive compounds, utilizing the 90 

26 acetylcholinesterase inhibitor (Lenz et al, 2012). The efficacy of Uncited references Q7 9 ^ 

27 AChEIs to reverse slow oscillatory activity in EEG has been 92 

28 demonstrated in AD patients. A consistent reduction of EEG delta Atrj et al 20Q8 . Babiloni et al _ 20 04; Babiloni et al., 2013 ; Birks, 93 

29 and theta oscillations and shift towards normalization was found 2006; Buckner et al 2 005; Carlen et al., 2012; De Graaf and 94 

30 after short as well as long term treatment with rivastigmine, Bonnard| 20 n ; Doody et al., 2001; Gray et al., 1989; Huxter et al., 95 

31 donepezil and tacrine (Adler and Brassen, 2001 ; Adler et al., 2004; 2003; Hyman et al _ 1990; , versen ig97; Jacobs et a , 2007; Jones 96 

32 Balkan etal., 2003; Brassen and Adler, 2003; Kogan et al., 2001). In and wilson 2005; Korotkova et al > 2010; L ee et al., 2005; Li 97 

33 addition, convergent imaging data demonstrated hyper meta- et a , 200g . Mangialasche et al., 2010; Maurer and McNaughton, 98 

34 bolism activation of brain structures known for their implication 20QJ . 0 ,, <eefe and Recce 19g3; prichep et a| 2006; Principles of 99 

35 in memory and cognitive functions such as fronto-parietal and Neural Sciencei 200 0; Reeves et al., 2002; Roberson and Mucke, 100 

36 hippocampal structures (Bohnen et al., 2005; Stefanova et al., 2006; Tallon _ Bau dry, 2009; Tallon-Baudry and Bertrand, 1999; 101 

37 2006). Uhlhaas and Singer, 2006; Weingartner, 1985. 102 

38 In the present work, all cognition enhancers consistently 103 

39 attenuated scopolamine-induced peak coherent delta activity in 104 

40 fronto-parietal networks. In line with previous results in AD pa- Acknowledgments 105 

41 tients showing the efficacy of AChEIs to attenuate the aberrant EEG 106 

42 network shift towards slow oscillations, the shift of EEG oscillations The authors are verv S rateful t0 Dr - J ason Wadde " and Dr - Tom 107 

43 towards high frequency oscillations is presumably suggestive of Van De Casteele for their statistical support. 108 

44 network neuroprotective activities for cognitive ameliorating effect 109 

45 in the scopolamine model. References 110 

46 111 

47 4.6. Slow theta (4-6.5 Hz) and gamma (32-48 Hz) oscillations and Adler ' G - Brassen, S., 2001. Short-term rivastigmine treatment reduces EEG slow- 112 
. „ ...... r . . wave power in Alzheimer patients. Neuropsychobiology 43, 273—276. . . „ 

48 connectivity: implications for mnemonic processing Ad , er G Brassen s chwalek, K., Dieter. B., Teufel, M., 2004. Prediction of treatment ! 1 3 

49 response to rivastigmine in Alzheimer's dementia. J. Neurol. Neurosurg. Psy- 114 

50 Mnemonic processing requires a cross-interaction between chiatry 75, 292-294. 115 
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51 hippocampal and cortical areas, while increased theta and gamma effectiveness of combination therapy in Alzheimer' disease. Alzheimer Dis. 116 
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55 be a fundamental mechanism by which the brain dynamically Zanetti, o., Zappasodi, E, Rossini, p.m., 2004. Mapping distributed sources of 120 
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